Theory of Brownian motion of Massive Particle in a Space with Curvature and Torsion 

and Crystals with Defects 



in 

On 
On 



H. Kleinert and S.V. Shabanov* 
Institut fur Theoretische Physik, Freie Universitat Berlin, Arnimallee 14, D-14195 Berlin 

(February 1, 2008) 

We develop a theory of Brownian motion of a massive particle, including the effects of inertia 
(Kramers' problem), in spaces with curvature and torsion. This is done by invoking the recently 
discovered generalized equivalence principle, according to which the equations of motion of a point 
particle in such spaces can be obtained from the Newton equation in euclidean space by means of 
a nonholonomic mapping. By this principle, the known Langevin equation in euclidean space goes 
over into the correct Langevin equation in the Cartan space. This, in turn, serves to derive the 
Kubo and Fokker-Planck equations satisfied by the particle distribution as a function of time in 
such a space. 

The theory can be applied to classical diffusion processes in crystals with defects. 



1. The classical equation of a massive point particle in a thermal environment reads 

^: ,,,:,■■', = /; + /;. ui 

o 

where p t is an arbitrary time-dependent external force and j\ a stochastic force caused by the thermal fluctuations 
(we use subscripts for the time variable) . The stochastic force may be modeled by a bath of harmonic oscillators of all 
frequencies u> as ]\ = J °° duX^X^. The oscillator coordinates satisfy the equations of motion X^ t + uP^X 1 ^ — A W ±J, 

y—{ • the right-hand side arising from the back-reaction of the particle. Solving the latter equation with respect to X^ t we 
! find ft as a functional of x\. Assuming an equal coupling of all oscillators, X u = y^r/Jn, we obtain the stochastic 

^""T 1 ■ differential equation 

o 
in 



mxi+ 7 xi-fi=4, (2) 
where rf t is called the noise variable. This has the decomposition 



2h ' 00 

S \ rf t — J doj\u X l u cos cut - sin cut . (3) 

g; 

For any given phase space distribution p B = p B (X,X) of initial oscillator velocities X^ and coordinates X^, Eq.(||) 
becomes a classical Langevin equation with noise averages being defined as mean- values with respect to the distribution 
J> . p B . In thermal equilibrium, p B follows the Boltzmann law ~ exp(— H B /kT), where the bath Hamiltonian is a sum 
of the oscillator energies H B = J duj{±l + uj 2 XI)/2. In this case, the noise is Gaussian and completely specified by 
its vanishing expectation and its two-point correlation function: 



(Vl)=0, (44) =67*r^'<y(*-0 • ( 4 ) 

2. By a nonholonomic mapping dx l — e t fl (q)dq^ jjj, the euclidean space is carried into an arbitrary metric-affinc 

The geometry in this space is defined by the metric — e^,e l „ and the 
The torsion is carried by the antisymmetric part of the affine connection, 
S^ x = (T^J" — T l/fl x )/2, which is a tensor, in contrast to itself. The curvature tensor arises from noncommuting 
derivatives of e l y. R^ V \ K — ei K (d^d l/ ~ df Jl d u )e t \. Thus, in a curved space, the nonholonomic quantities ei^(q) fail to 
satisfy Schwarz' criterion and are no proper functions p|. 



space with curvature and torsion 
affine connection V ^ v x — e i x d lJL e' l l 
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Recently it was found that the nonholonomic mapping dx l — e l ^(q)dq fl transforms euclidean equations of motion, 
Schroedinger equations, and path integrals correctly into spaces with curvature and torsion This universal mapping 
rule was called quantum equivalence principle (QEP)H. 

Applying this principle to the equation of motion (|2|), we obtain the analogous equation in the general metric- affinc 
space: 

m (qt + T^qUt) + 7# - ft = ^{q t )4 ■ (5) 

To obtain physical consequences, we must find equations in which the nonholonomic quantities e^^qt) are eliminated 
in favor of the well-defined geometrical objects, the metric g^iqt) and the affine connection r Atl/ A (q , t ). This is possible 
by deriving from Kubo's stochastic Liouville equation and the Fokker- Planck equation. 

For this purpose we rewrite the Langevin equation as a system of two first-order differential equations 

# = -g^ivtH ; (6) 

m 

= -Fl + ^\. (7) 

The last equation defines the total apparent force F*, where /* = gfj,v(qt)ft- We mark the time dependence also 
by a superscript to avoid a pile-up of subscripts. Note that although the time derivative of the momentum p 1 ^ and 
the force do not transform as vectors under general coordinate transformations, Eq. ^ is nevertheless covariant 
under these transformations, just as the Langevin equation (|5|). 

At each time £, the system following Eqs. (^) and (^) is in a microscopic state with the distribution function 
S(p — p t )S(q — qt) (omitting the spatial indices p). This can be thought of as a conditional distribution function 
determining the distribution of p and q for a given initial distribution function 5{p — p°)5(q — qo)- If the initial values 
of p l and qt are distributed with the probability density p(p a , qo), the distribution function at any later time t can be 
found by an average over the initial distribution 

p?(p, q) = [ dp°d q() p(p , q )S(p ~ p*).5(g - qt) , (8) 



where qt = qt(p°, qo) and p l = P t (p°, qo) are solutions of the system (^|), Q with the initial conditions p* _0 = p° and 
qt=o = qo- Upon taking a time derivative of (^) and making use of both (g) and (Q) and of the identity 

U{z - z t ) = z®S(z - z t ) = [z t S(z - z t )} = ~ [V(z)S(z - z t )} , (9) 
dt ozt oz oz 

where a dynamical variable Zt (any of p^ or q^) obeys the equation of motion it = V(zt), we find Kubo's stochastic 
Liouville equation 

ftp? = L(rk)p? , (10) 
with L(rjt) being the noise-dependent Liouville operator 

Hm) = ~o^q)^ + ±.o[F ll (p tq )-e* l rS\ . (11) 

The symbol o emphasizes the fact that the differential operators in front of it act on the product of the function 
behind it with p\ in Eq. @. 

It is straightforward to verify the invariance of the Liouville operator and, hence, of Kubo's equation with respect 
to general coordinate transformations. 

A solution of Kubo's stochastic equation (|o|) is a noise-dependent distribution function which determines the 
probability to find a particle in an infinitesimal volume dq by 

dPt(q) = dqj dp(pUp, q)) = dq J dpp t (p, q) . (12) 

It follows then from (||) that J dPt(q) = J dPo(q) = 1, i.e. the temporal evolution of the probability distribution 
described by Kubo's equation (|o|) is unitary. 
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The locality of the noise correlator ([|) enables us to derive a Fokker-Planck equation governing the temporal 
evolution of the noise-averaged distribution p t (p, q). For this let us first calculate the average ip\ — {if t 5{p—p t )8{q— q t )) . 



Generalizing the partial integration formula J e ar) l 2 r\]{r])dr] = a 1 J e aT) ^ 2 f'(rj) = —a 1 J 
Gaussian noise, we find 



f{rj) to any 



8rf t , 



5(p-p t )6(q - q t ) 



(13) 



where we have used the explicit form of the two-point noise correlation function (Q) and the identity (0). To calculate 
the variational derivatives of dynamical variables with respect to noise in (|l3|), we formally integrate (p|) and (J7]): 



Pl=pl + J dt' (-i^' + e^,) 



q? = q% + - / dt'g^(q v )pl 



(14) 
(15) 



Taking the variational derivative of (|l4|), we obtain 



dt" 



SF, 



t" 



8r& 



dt"e?(q t „)S(t"-t') 



(16) 



The Langevin equation (Q) is causal which implies that p^,, or q^,, depend on rf t , only for t" > t' . This has been used 
to restrict the integration range in (|l6|). The first integral in the left-hand side of (|lj) tends to zero as t' approaches 
t, whereas the second integral is equal to ef(qt')Q(t — t') and, hence, is not uniquely determined, since it involves the 
value of Heaviside function at zero argument, 8(0), when t = t' . Therefore a regularization is needed. We replace 
the delta-function in the correlator (||) by a smooth would-be delta- function of width e, S e (t — t'), which satisfies 
dt6 e (t) = 1 and <5 e (£) = 5 e (— t). It is easily seen from (||) that such a regularization can be achieved by retaining 



a weak dependence on oj in the coupling constant 
by 



With this regularization, we have to replace Sp l /Sr]l in (Hi 



dt'S £ (t 



- = J dt'S e {t 



t')e l (q t ,)0(t-ti) = -e*( qt ) , 



(17) 



where we have dropped the contribution of the first integral of ( |iq ) since it vanishes for e — > 0. 

Considering analogously 8q^ /5rf t , : we conclude that the latter variational derivative vanishes as t' approaches t. 
Averaging tf\ with respect to the initial distribution p(po, <?o) we nn d the following relation 



f r 9 

j dp dq p(p ,qo)(pl = {rfj.pl) = -37*;Te^- — p t 



(18) 



Taking the noise average of Kubo's stochastic equation (|To|) and making use of (|l^), we end up with the Fokker-Planck 
equation associated with the Langevin equation (^|): 



d t pt = L T pt , 
Lt — 



- o 

p m 



dpp 



F ll (p,q)+3 1 kTg^ 



(19) 
(20) 



1 This physical regularization corresponds to the Stratonovich treatment of stochastic integrals 
smoothed correlator appears also due to quantum effects (the fluctuation-dissipation theorem) . 



p. 55. A symmetrically 
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Thus, we have eliminated the nonholonomic mapping functions in favor of the well-defined affine connection and 
metric. Integrating (19) over all phase space we arrive at the probability conservation law: d/dt J dpdqpt = 0, where 
we have used the fact of vanishing surface integrals occurring upon the integration of the right-hand side of ([jl]) . 

In the overdamped limit, we may drop the mass term in (^|) and derive analogously an equation for the g-space 
probability distribution D t (q), defined by dPt(q) = dq^JgD t {q), 



d t D t = 



3kT 



<9 M o ef9„ oe^y/g- 



3fcT 



d» o g^y/gidu + 2S V - nf v ) 



D, 



(21) 



where k = (3kTj) , S v = and g = detg^. Since S v and f v are vectors, the operator in the brackets is 

invariant under general coordinate transformations. In the absence of torsion, the limiting equation (^lj) agrees with 
a result in Ref . fjn . 

3. Since defects in crystals can be described geometrically by a nonvanishing torsion and curvature in a metric- affine 
space in which distances are measured by counting atomic steps |jj , the Fokker-Planck equation (^) can be applied 
to the classical diffusion of interstitial atoms in crystals M . 

It must be realized, however, that the approximation is necessarily very crude. The reason for this is the twofold 
role played by the phonon bath. First, it provides the system with a random driving force which in our equations 
has been simulated by the model noise r\\. Second, it influences the hopping probabilities of the interstitials. The 
latter effect has not been accounted for in our initial Eq. (||). In crystals, the hopping probability depends on 
elastic distortions which correspond to time-independent general coordinate transformations in the above theory. The 
Langevin equation (||) is, however, invariant with respect to these transformations. A proper treatment should take 
into account the fact that elastic distortions are also produced by the off-shell phonons comprising the noise. 

Another field of applications should be positrons or muons in a crystal at low temperatures below a few or about 
80 Kelvin, respectively. These particles are bound in the lowest atomic orbits, so that their wave functions follow 
the tight-binding approximation. Spending most of their time near the atomic sites, they live effectively in the above 
crystal geometry. 



It is worth noting that in the absence of external forces, the overdamped equation (21) agrees with an equation 
recently studied in Ref. H. 
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